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Abstract: We study contributions of a scalar diquark particle in a color representation 
of anti-triplet and sextet to the top quark pair production at the Tevatron and the Large 
hadron collider (LHC). The model can give Forward-Backward (FB) asymmetry at the 
Tevatron while can avoid the same sign top quark production at the LHC by assuming the 
top-number conserving diquark couplings. We study compatibility between the large posi- 
tive FB asymmetry observed at the Tevatron and non-observation of the charge asymmetry 
at the LHC, by including contributions from the single and pair production of diquarks. 
We find that the whole parameter space of the models can soon be explored at the LHC by 
measuring the total ti production cross section and the inclusive charge asymmetry with 
smaller uncertainties. In addition, we compare the statistical significance of the charge 
asymmetry measured at the LHC with that of the optimal observable of the sub-process 
FB asymmetry, and find that they are comparable even when we ignore the uncertainty in 
the parton distribution functions. 
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1 Introduction 

The charge asymmetry in the ti production is the difference between the top and anti-top 
quark distributions. In the standard model (SM), small charge asymmetry is expected in 
the light quark pair annihilation processes, qq — > ti at the next-to-leading order (NLO) of 
QCD perturbation theory [1], which predicts small Forward-Backward (FB) asymmetry in 
pp collisions (Tevatron) and tiny charge asymmetry in pp collisions (LHC). 

At the Tevatron ti pair is produced mainly in the quark pair annihilation process, 
qq — > ti, where the quarks (anti-quarks) are mainly moving along the proton (anti-proton) 
momentum direction. The FB asymmetry can be defined as 



N(Ay > 0) - N(Ay < 0) 
FB N (Ay > 0) + N(Ay < 0) ' { ' 

where N is the number of events and Ay = yt — Ut is the difference in rapidity of top and 
anti-top quarks along the proton momentum direction in the laboratory frame. 

Recent results from the CDF and the DO collaboration at the Tevatron report positive 
asymmetries [2, 3] 

A^ F = 15.8 ± 7.4 % (1.2a) 
= 19.6 ± 6.5 % (1.2b) 

at the sub-process level after correcting for backgrounds and detector effects, while the 
SM prediction at the NLO in QCD is 5.0 ±0.1 % [4-6]. Furthermore, the CDF reported 
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even larger asymmetry of 47.5 ± 11.4% at A% > 450GeV [2] while the DO reported the 
asymmetry of 15.2 ± 4.0% in the rapidity of leptons [3], both exceeding the SM predictions 
by more than 3 standard deviations. 

At the LHC, the sub-process FB asymmetry cannot be observed from Ay = yt — Vt 
distribution because it is app collider. Instead, the following charge asymmetry is sensitive 
to the qq — > ti sub-process FB asymmetry at the LHC 

= N(A\y\>0)-N(A\y\<0) 
C N(A\y\ >0) + A(A|y| < 0) ' 1 ' 

where A\y\ = \yt\ — is the difference in the absolute values of rapidities of top and 
anti-top quarks. The latest results from the ATLAS and the CMS collaboration at the 
LHC are [7, 8] 

^atlas _ _ L8 ± 2 g ^ stat j ± 2 3 ( sygt _) %) ^ 4a ^ 

A CMS = _ 13 ± 2 8 ( stat ) +2.9 (sygt ) %) (L4b) 

which are within the uncertainty consistent with the SM prediction, Ac = 1.15±0.06 % [9]. 
The asymmetry is small because the gluon fusion sub-process, gg — > ti, gives the dominant 
contribution at the LHC. Nevertheless we expect the LHC data to reveal the asymmetry 
with the help of its high energy and its expected high luminosity. 

In this paper, we study contributions of a scalar diquark particle in a color repre- 
sentation of anti-triplet and sextet to the total ti cross section and inclusive asymmetries 
at the Tevatron and the LHC. We evaluate the parameter space of the models which are 
consistent with the large positive FB asymmetry at the Tevatron, non-observation of the 
charge asymmetry at the LHC and non-deviation of the total ti production cross section 
from the SM prediction both at the Tevatron and at the LHC. 

The contribution of diquarks to the FB asymmetry at the Tevatron is discussed in 
refs. [10-23] and their contribution to the charge asymmetry at the LHC is discussed in 
refs.[24, 25]. The (j) single and pair production expected in the models is discussed in 
refs. [10-12, 14, 15, 17, 19-22, 24, 26]. We explore the parameter space of the models 
by including contributions from the single and pair production of diquarks to the charge 
asymmetry Aq at the LHC. 

This paper organized as follows. In the next section, we introduce the diquark models 
and discuss their phenomenological consequences at the Tevatron and the LHC. In Section 
3, we show our numerical results of the diquark models. In Section 4, we introduce the 
optimal observable of the sub-process FB asymmetry. In Section 5, we compare the sta- 
tistical significance of the charge asymmetry Aq eq.(1.3) measured at the LHC with that 
of the optimal observable of the sub-process FB asymmetry, and examine the efficiency of 
Aq- The last section gives the summary. 
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2 Diquark model 



We consider a model which consists of a new scalar boson </> in a SU(3) color representation 
of anti-triplet or sextet with a diquark quantum number: 

C = \V2C^t R ■ d j R 4> a * + h.c, (2.1) 

where C?- is Clebsch-Gordon coefficients and described clearly in the appendix of ref.[27]. 
The representations under SU(3)c x SU(2)l x U(1)y for anti-triplet and sextet diquarks 
are (3,1,1/3) and (6,1,1/3), respectively. We may denote the anti-triplet diquark as 
03 = (cj) 1 , cj) 2 , (f> i ) T and the sextet diquark as (J)q = (cj) 4 , (j) 5 , <p G , (j) 7 , <p s , <ft 9 ) T ■ If only the 
interactions of eq.(2.1) are present, the diquarks carry the quantum number of the top 
and down quarks and their chiralities, hence each flavor number as well as the chirality is 
conserved. 

The SM gauge invariance allows top-light diquark couplings between the doublets 
(^Li^l) and (uL,di), and between the singlets t R and d R , or u R . The diquarks of t R and 
u R have the same form of the Lagrangian eq.(2.1) when d J R is replaced by u° R and <j) a have 
the hypercharge 4/3. The SU(2)l representation of the {tL,bi) — (ul,cIl) diquarks is either 
singlet or triplet, (j) or <fj' = (<j>" 4 / s , 4>"i/ 3 , ^'_2/3) T ' respectively, whose Lagrangians are 

£su(2) L singlet = Vv^Cg-^^ • d{ - b\ ■ u\)^* + h.c, (2.2a) 

^SU(2) L triplet = A" V2C°j l-t L ■ U{4>1% + b{ ■ dy_ a * /3 + i=(4 • d{ + b\ ■ 4)^*1 + h.C. 

(2.2b) 

The electromagnetic (EM) charge 1/3 diquarks contribute to dd —> tt, whereas the EM 
charge 4/3 diquarks contribute to uu — > ti. We find that all the diquarks give positive FB 
asymmetry, especially at high M t j , because of the chirality conservation of each Lagrangian. 
The left-chirality diquarks of eq.(2.2) also contribute to the single top quark production, 
ud — > tb, and hence are strongly constrained [28, 29]. 

Rather than studying all of them exhaustively, we choose to study the t R -d R diquarks 
of eq.(2.1), since the anti-triplet diquark has the least exotic quantum number among 
all the diquarks, which can be the super-partner of the right-handed down quark if the 
supersymmetry breaking scale is far above the electroweak scale [30] . 

The diquark interactions of eq.(2.1) can generate the sub-process FB asymmetry at 
the leading order in the process of 

diipi, Ai) + dj(p 2 , A 2 ) ->■ t k (p 3 , A 3 ) + ii(p 4 , A 4 ) (2.3) 

through the it-channel exchange of 4>, as illustrated in Figure 1, where pi and Aj are momenta 
and helicities, respectively, and i, j, k and / are the color indices. 
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Figure 1. Feynman diagrams of the process dd — > ti through the u-channel exchange of 0| 6 , in 
addition to the QCD amplitude with s-channel exchange of gluon (g a ). Hence i, j, k, I and a are 
the color indices, p n and X n are momenta and hclicities, respectively. 



The helicity amplitudes are summarized below, 

= E CuC^rm sin + ^ £ T^l^rrH sin 9 (2.4a) 



a a=l 



<t- = ^ E C ^ (4) 2 (1 + W + cos e ) + 4 -l± ipt (f) 2 (1 + cos 0) 



9 a \ / a=l 



(2.4b) 



cos I 



(2.4c) 

< + v = -^E^^ si ^-^E^f^ sin0 > ( 2 - 4d ) 

a o=l 

where the first terms correspond to the amplitudes of the u-channel exchange of 0| or <^g 
and the second terms give the QCD one-gluon exchange amplitudes. Initial down quarks 
are assumed massless, and hence the diquarks do not contribute to the M^ 3 ^ 4 amplitudes, 
which are obtained from the QCD part of the amplitudes eq.(2.4) by parity transformation 
(Aj — > — Aj, cos 9 — > — cos 9). Here s = (j>\ + P2) 2 is the invariant mass squared of the ti 
system, f3 = (1 — 4m|/s) 1//2 is the velocity of the top quark, 9 is the polar angle between 
the initial down quark and the final top quark momenta in the ti rest frame. 

The color-space propagators of the 3 and 6 diquarks as well as that of QCD gluons 
are, respectively, 

3 a a* 1 

E C il C jk = 2^ji 5 kl ~ SjlSki), ( 2 - 5a ) 
o=l 

9 1 

E c ii c Tk = 2^ iSki + (2 - 5b) 

E T ji T kl = 2^ lSki ~ ( 2 - 5c ) 
a=l 
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and the color factors which appears in the color summed squared amplitudes are 



3 3 8 

E EEW^ = - 2 . ( 2 - 6a ) 

i,j,k,l=l o=l 6=1 
3 9 8 

E EEWWi = +2, (2.6b) 

i,j,k,l=l a=A 6=1 
3 3 

E E( c ^)*^- = 3 ' ( 2 - 6c ) 

i,j,k,l=l a, 6=1 
3 9 

E E( c ^ ! )* c ^ = 6 ' ( 2 - 6d ) 

i,j,k,l=l a,b=4 
3 8 

E E( T ^)* T ^ = 2 - ( 2 - 6e ) 

i,j,k,l=l a, 6=1 

At this stage, we can tell that the leading part of the QCD amplitudes interfere construc- 
tively (destructively) with the sextet (anti-triplet) diquark exchange amplitudes, and we 
expect positive FB asymmetry for the sextet contribution at all energies. Even for the 
anti-triplet, we find positive FB asymmetries for strong couplings A when the diquark ex- 
change amplitude dominates over the QCD amplitude in eq.(2.4b), which happens e.g. at 
= 0.46,0.33,0.14 for A = 2.6,2.8,3.0 with = 500GeV, respectively. The positive 
FB asymmetry is a consequence of the chirality conservation of the effective Lagrangian 
eq.(2.1), and hence it is common for all the diquark models that respect the SM gauge 
invariance as in eq.(2.2). 

At the LHC, the positive FB asymmetry at the sub-process level gives rise to the 
positive charge asymmetry Ac eq.(1.3). In addition, single and pair production of diquarks 
can be significant at the LHC, 

dg -)• 4>i ->• tid, (2.7a) 
dg^ i0 -> tid, (2.7b) 
gg — > <j)(j) —> tidd, (2-7c) 
qq ->• (j>4> ->• tidd. (2.7d) 



If the diquark <fi is heavier than the top quark, > m t , the single production processes 
give tt+jet events and the pair production processes give tt + 2jets events. These additional 
processes can also contribute to the inclusive FB and charge asymmetry as well as to the 
total tt production rate. The relevant diagrams for the single and pair production of 
diquarks via dg — > (j)t and dd — > (fxfi are shown in Figure 2 and Figure 3, respectively. 

The single <fi production process of eq.(2.7a) is dominated by the ^-channel top quark 
exchange amplitude in Figure 2(b) when the diquark is heavier than the top quark, and 
hence the diquark is produced mainly along the initial down quark momentum direction. 
This leads to yt > yi at the Tevatron. At the LHC, initial down quarks in protons are 
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d ( a ) d d (b) d 



Figure 3. Feynman diagrams of the process dd — > (jxfi. 

more boosted on average than gluons, and hence the diquark is mainly produced at large 
1 7/0 1, which leads to \y t \ > \yi\. 

Similarly in the single (j) production process of eq.(2.7b), the anti-diquark <j) is pro- 
duced mainly along the initial anti-down quark momentum direction due to the top quark 
exchange amplitude. This also leads to yt > yi at the Tevatron. At the LHC, however, 
since anti-down quarks have only small energy fractions of incoming protons, the charge 
asymmetry from the single <f> production is negligibly small. 

In the (j) pair production process of eq.(2.7d), the t-channel top quark exchange ampli- 
tude in Figure 3(a) gives a sub-process FB asymmetry y^ > y^, which leads to yt > yt at 
the Tevatron and \yt\ > \yt\ at the LHC. The qualitative effect of the charge asymmetry 
at the LHC is, however, negligibly small because of the smallness of the d energy fraction 
in the proton. The <j) pair production have contributions mainly to the total ti production 
cross section, mainly via the gg — > <j)(f) sub-process of eq.(2.7c) which gives no asymmetry. 

From the above discussion, single and pair production of diquarks also contribute 
positively to the FB asymmetry at the Tevatron and to the charge asymmetry at the LHC, 
as well as to the total ti production rate. 

3 Numerical results of the diquark models 

3.1 Total ti cross section and asymmetries at Tevatron and LHC 

In this section, we examine the parameter space of the diquark models and look for allowed 
regions which are consistent with the measurements both at the Tevatron and the LHC. 
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At the Tevatron, while producing a large positive FB asymmetry, it is important to 
keep the deviation of the total cross section of the tt production a(pp —> tiX) from the SM 
prediction small because it has been measured to agree with the SM prediction [31-35]. 
Thus we require the following two conditions to be satisfied, 

A FB > 5%, (3.1a) 
\<tsm+np-(t sm \ ^ < Q ^ 
&SM 

in the leading QCD order, and consider that the parameter region is excluded if one of the 
above conditions is not satisfied. Since the SM predicts Afb of about 5% [4-6] at the NLO 
of perturbative QCD, our requirement of Afb > 5% at the leading order may be large 
enough to be consistent with the large positive FB asymmetry observed at the Tevatron. 
The condition eq.(3.1b) for the total cross section is rather loose, since the observed tt 
production cross sections [31-35] tend to be somewhat larger than the NLO+NLL QCD 
prediction of about 7.14 pb for mt = 172.5 GeV [36]. 

At the LHC, on the other hand, since the uncertainties on the tt total cross section 
and the charge asymmetry Ac are still large, we use the observed upper bounds from the 
experimental data. As for the total cross section, we require [37] 

a(pp -»• ttX) < 204 pb at 95% confidence level (C.L.), (3.2) 

and for the charge asymmetry, we require one of the following constraints [7] 

A c < 2.8 % at 90% C.L., (3.3a) 
A c < 4.1 % at 95% C.L.. (3.3b) 

We use CTEQ6L1 [38] parton distribution function (PDF), and take the SM param- 
eters as the top quark mass m t = 172.5 GeV and the QCD coupling as(mz) = 0.130 
which is used in CTEQ6L1, throughout this paper. All the results are obtained by using 
MadGraph [39-41] at the matrix element level. We do not take into account the effects 
from parton showering, hadronizations and detector conditions. 

The result is shown in Figure 4. The horizontal axis corresponds to the diquark mass 
and the vertical axis to the coupling constant A. The anti-triplet diquark model is 
evaluated in the top panels, and the sextet model is in the down panels. The two panels on 
the left hand side are obtained without contributions from the single and pair production 
of diquarks, and the two figures on the right hand side are obtained with contributions 
from the single and pair production of diquarks. Two dashed lines reflect the conditions 
at the Tevatron eq.(3.1), and the shaded region satisfies both of these two conditions. The 
space on the right hand side (RHS) of the solid line satisfies both eqs.(3.2) and (3.3b), and 
the space on the RHS of the dotted line satisfies both eq.(3.2) and the stronger constraint 
eq.(3.3a). 

Comparing the left and the right panels in Figure 4, we find that the contribution from 
the single and pair production of diquarks is quite significant. The allowed region of the 
anti-triplet model is significantly reduced when (p production contributions are taken into 
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Anti-triplet w/o <p productions Anti-triplet with (j> productions 




Figure 4. The horizontal axis corresponds to the diquark mass and the vertical axis to 
the coupling constant A. The anti-triplet diquark model is evaluated in the top panels, and the 
sextet model is in the down panels. The two panels on the left hand side are obtained without 
contributions from the single and pair production of diquarks, and the two panels on the right hand 
side are obtained with contributions from the single and pair production of diquarks. Two dashed 
lines reflect the conditions at the Tevatron eq.(3.1), and the shaded region satisfies both of these 
two conditions. The space on the right hand side (RHS) of the solid line satisfies both eqs.(3.2) 
and (3.3b), and the space on the RHS of the dotted line satisfies both eq.(3.2) and the stronger 
constraint eq.(3.3a). 

account in the right panel. It is also worth noting that </>•) production contributes to the 
positive FB asymmetry even at the Tevatron when < 400 GeV. To a lesser extent, 4>q 
production contributions reduce the allowed region of the sextet model as shown by the 
bottom two panels in Figure 4. We should therefore consider contributions from the single 
and pair production of diquarks, as well as from the u-channel diquark exchange process, 
when we discuss the FB and charge asymmetry in the diquark models. 

The anti-triplet diquark model can be excluded if we adopt the constraint on Ac at 90% 
C.L. from the LHC eq.(3.3a), whereas the sextet diquark model still has small compatible 
parameter space in the RHS of the dotted line in the shaded region of the bottom-right 
panel of Figure 4. If we loosen the constraint on Ac to 95% C.L. limit of eq.(3.3b), both 
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models still have allowed parameter spaces. Nevertheless, the allowed parameter regions 
are rather narrow, and can soon be explored by measuring the total ti production cross 
section and the inclusive charge asymmetry with smaller uncertainties at the LHC. 

3.2 Acceptances of the ti events induced by diquark productions 

In Section 3.1, we find that effects of the single and pair production rate of diquarks can be 
large at the LHC. When we discuss these contributions to Ac and a(pp — > tiX) in Figure 
4, we assume that 100% of the diquark production contribute to the inclusive ti events. 
However, in actual measurements of Ac and cr(pp —> tiX), experimentalists apply event 
selection cuts on final states in order to increase the ratio of the number of the signal ti 
events over that of background events. Since the single and pair productions of the diquark 
lead to ti+ 1 jet and ti+2 jets events, respectively, at the sub-process level, the probability 
of those events to pass the selection cuts can differ from that of the SM it events. 

In this section, we study the acceptance of the ti events induced by the single and pair 
productions of the diquark by using Monte Carlo generated event samples. For this study, 
MadGraph/MadEvent [39-41] is used to obtain the parton level distributions for ti+n 
jets events, which are interfaced to Pythia6 [42] for the parton showering with the shower 
kT matching scheme [43]. We generate ti events where one of the W bosons decays into 
an electron or muon and the corresponding neutrino and the other W boson decays into a 
pair of jets. After jet clustering, we require jets to satisfy, 



p J T > 30 GeV, 
| rf\ < 2.5, 
Number of jets > 4. 



(3.4a) 
(3.4b) 
(3.4c) 



We further impose the following selection cuts, 



EJp iss > 35 GeV, 
m T (l, v) > 25 GeV, 



p e T > 30 GeV, 
7] e \ < 2.5, 



(3.5a) 
(3.5b) 
(3.5c) 
(3.5d) 



-Total 
Jet 




(3.5e) 



for the electron channel, and 



p T > 20 GeV, 
if \ < 2.1, 



(3.6a) 
(3.6b) 
(3.6c) 
(3.6d) 



Ef lss > 20 GeV, 
Ef iss + m T (l,u) > 60 GeV, 




Pt 



(3.6e) 
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^^_rn [GeV] 
Process 


200 


400 


600 


800 


1000 


(j)t, t(f) 


1.05 


1.11 


1.15 


1.15 


1.14 


# 


1.07 


1.17 


1.18 


1.18 


1.18 



Table 1. Acceptances of the single and pair productions of the diquark normalized to that of the 
SM tt process for different masses of the diquark. 



for the muon channel, where tut (I, v) is the transverse mass J2p l T E™ iss (l - cos(0' - ^ miss )) 

and E]°^ l {AR < 0.4) is the sum of jet energies in a cone with AR = y/Atfj + Ar/2. = 0.4 
around the lepton track. These event selection cuts are obtained from refs.[7, 8], although 
we do not simulate detector effects in this paper. We simply define the acceptance as a 
ratio of the number of events that pass the above selection cuts eqs.(3.4, 3.5, 3.6) over the 
number of generated events. The acceptance for the SM ti events are found to be 0.45 in 
our analysis. 



We show in Table 1 the acceptances of the single and pair production processes nor- 
malized to that of the SM ti process for different masses of the diquark at y/s = 7 TeV. As 
the mass of the diquark grows, we expect that not only the down-quark jet but also the top 
quark from the diquark decays are boosted, resulting in hard jets and a hard lepton in the 
final state which make it easy for the event to pass the selection cuts. However, we find in 
Table 1 that the acceptance does not simply increase as the mass grows. The reason is as 
follows. For heavy diquark with mass > 600 GeV, the top quark from its decay is highly 
boosted and when this top quark decays leptonically, a lepton and a bottom quark tend to 
be collinear. These events tend to be rejected by the present selection cuts for the lepton 
isolation, eqs.(3.5e, 3.6e). We have confirmed that the acceptances increase almost linearly 
with the diquark mass when the lepton isolation cuts, eqs.(3.5e, 3.6e), are removed. All the 
results in Table 1 are obtained for the color anti-triplet diquark. We find similar results 
for the color sextet diquark. 



In order to take into account higher acceptances for the diquark production events in 
the analysis of Figure 4, we simply multiply the diquark production cross sections by the 
acceptance ratios in Table 1. The cross sections become larger and therefore the solid and 
dotted lines shift a bit toward the right in the two panels on the right hand side in Figure 
4, resulting in slightly reduced allowed regions. However, because the acceptance ratios 
are only 5 to 20% larger than the unity and also because the allowed regions in Figure 4 
corresponds to the regions where the diquark production cross sections are small, we find 
that the results presented in Section 3.1 are not affected significantly by taking account of 
the difference in the acceptance of the diquark production events and the SM ti events. 
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4 Projection method 



4.1 General formula 

In this section, we would like to study the relationship between the sub-process FB asym- 
metry and the charge asymmetry Ac, eq.(1.3), which is measured at the LHC. 

The number of top and anti-top pair production events in pp collisions can be described 

as 

N( T ,Y,y) =c\D gg {T,Y)a gg (T,y) + D q q(T,Y)a q g(T,y) + Dq q (T,Y)aq q (T,y)} (41a) 



X 



D gg (r, Y)a gg (r, y) + Df(r, Y)af(r, y) + £>«(r, Y)af{r, y) , (4.1b) 



where C is an integrated luminosity, a a b(T,y) are the ab — > ti sub-process cross sections, 
r = X1X2 gives the product of the momentum fractions of the colliding partons, y is the 
rapidity of the top quark in the sub-process rest frame, and Y is the rapidity of the ti 
system in the laboratory frame. In eq.(4.1a), D ab {r,Y) are the products of PDFs, 

D gg (r, Y) = D g/p (x a ) x D g/p (x b ) = D g/p {^e +Y ) x D g/p (^/re~ Y ), (4.2a) 
D qq (r, Y) = D q/p (x a ) x D q/p {x b ) = D q/p (VTe+ Y ) x D q/p {^ e - Y ), (4.2b) 
D qq (T,Y) = D q/p (x a ) x D q/p (x b ) = D q/p (^e +Y ) x D q/p (^ e - Y ). (4.2c) 

In eq.(4.1b), a Q g and a q ^ are, respectively, the symmetric and anti-symmetric part of the 
sub-process cross section about y, or equivalently symmetric and anti-symmetric under the 
exchange of colliding q and q momentum direction, 

<7qq(T, y) = Vqq(T, ~y) = of^T, y) + CT^(t, J/), (4.3a) 
a qq{ T , -V) = Vqqij, V) = <4 ? ( r > 2/) ~ °"a( T ' V)' ( 4 ' 3b ) 

Similarly, D q ^ and D q £ in eq.(4.1b) denote the symmetric and anti-symmetric combinations 
of Dgq(r,Y) in eq.(4.2b) and D qq (r,Y) in eq.(4.2c), 

Df(r, Y) = D qq (r, Y) + Dqq(r, Y), (4.4a) 
Df (r, F) = D M -(r, Y) - D qq (r, Y). (4.4b) 

In Figure 5, we show the Y distribution of the product of the d and d quark PDF's; the 
black solid curve gives D d j(r,Y) of eq.(4.2b), the black dashed curve gives D^ d (r,Y) of 
eq.(4.2c), whereas the blue and red curves are their symmetric Dg d (T,Y) of eq.(4.4a) and 
anti-symmetric D^{t,Y) of eq. (4.4b), respectively. 

In eq.(4.1b), multiplying both sides by D q ^{T,Y) and integrating about Y in its whole 
region, the first and second terms on the right hand side vanish because these terms become 
anti-symmetric about Y, and only the third term survives, 

f dY r +1 drD q y q (T,Y)N(r,Y,y) = £ [ dY dT(D qq (r,Y)f a qq {r,y). (4.5) 

J all J r, J all J r, 
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Figure 5. The Y distribution of the product of the d and d quark PDF's at 7 TeV LHC; the black 
solid curve gives D d( j(T, Y) of eq.(4.2b), the black dashed curve gives Z?j d (r, Y) of eq.(4.2c), whereas 
the blue and red curves are their symmetric Dg(r,Y) of eq.(4.4a) and anti-symmetric D a (t, Y) 
of eq.(4.4b), respectively. The product of the momentum fractions of the colliding partons, r, is set 
to t = 4m,(/s, and the factorization scale is set to m t . CTEQ6L1 PDF is used [38]. 

In other words, the asymmetric part of the sub-process cross section, a A (r,y), can be 
projected out even in pp collisions if we know the q and q PDF accurately. In the following 
sections, we examine the possibility of the optimal measurement of the asymmetric cross 
section of the qq — > ti sub-process in the absence of the PDF uncertainty. 

Making an approximation that cr q A {Tj) is constant between tj and Tj+i, we obtain the 
asymmetric part of a sub-process cross section a q ^{rj) as follows 



Since there is no statistical correlation between N(T,Y,y)6(y) and N(T,Y,y)9(—y), the 
statistical uncertainty of & A (rj) is estimated as, 



If we can ignore the PDF uncertainty in D a (t,Y), the projected cross section eq.(4.6) 
should have the smallest error since the weight D A (r,Y) minimizes its statistical uncer- 
tainty [44-47]. 

4.2 Extracting a A q (r) from MC generated event samples 

In this section, we demonstrate the extraction of ct a (t) from Monte Carlo (MC) generated 
ti event samples by using eq.(4.6). Event samples are generated by MadGraph/MadEvent 
[39-41] assuming the color anti-triplet and sextet diquark models. FeynRules is used to 
make the model files for MadGraph [48]. From the allowed parameter spaces in Figure 
4, we choose two points, (m^ 3 , A) = (500,2.65) and (600,2.95) for the anti-triplet model, 
(m^ 6 , A) = (700, 2.1) and (1000, 2.75) for the sextet model. Since the diquark model eq.(2.1) 
gives the sub-process FB asymmetry at the leading order in the process of dd — > ti eq.(2.3) 
through the it-channel exchange of the diquark, only the dd contribution, <t^(t), is non- 
zero at the leading QCD order, and this is what we try to extract from MC generated 




Sail dY % +1 dr f all dy [%) - 8{-y)] Pf{r, Y)N(r, Y, y) 
£ LldY S T 4 +1 dr(Df(r,Y)) 2 



(4.6) 




y/lau dY ST dr J all dy [%) + 9(-y)] (iff (r, Y)) 2 N(r, Y, y) 
£ LldY jj+i dr{Df{r,Y)f 



(4.7) 
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Figure 6. The results of extracting aj^ from MC generated ti event samples as a function of the 
ti invariant mass M t (, with the statistical uncertainties. The theoretical calculations arc shown as 
red solid curves. The color anti-triplet (top panels) and sextet diquark (bottom panels) models arc 
assumed. Parameter sets for the diquark mass m<f, and the coupling constant A are shown above 
each panel. 



ti event samples by using eq.(4.6). The MC generated event samples correspond to an 
integrated luminosity of 20fb~ 1 at the 7 TeV LHC and we take into account the overall 
QCD NLO K factor of 1.5 and the branching ratio 0.29 of the semi-leptonic ti decay into 
Z^+jets for I = e or fi. 

The results of extracting ct^(t) with the statistical uncertainties are shown in Figure 
6. The theoretical calculations of (r) are shown as red solid curves. The horizontal 
axis is the ti invariant mass, M f j = ^/TS. The results are consistent with the theoretical 
calculations within the statistical uncertainties, and hence we confirm that the projection 
formula eq.(4.6) works. 



5 Evaluation of the charge asymmetry Aq 

The charge asymmetry Aq of eq.(1.3) measured at the LHC uses the difference A|y| = 
| yt | — | yi | in the magnitudes of the t and i rapidities along the proton momentum direction 
in the laboratory frame. In terms of Y and y, \yt\ and \yi\ are written as 



\y t \ = \Y + y\ = {Y + y)9{Y + y)- 
\Vi\ = \Y -y\ = (Y -y)9(Y-y) 



(Y + y)8(-(Y + y)), 
(Y — y)6(—(Y — y)). 



(5.1a) 
(5.1b) 
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Therefore, the numerator of the charge asymmetry eq.(1.3) can be expressed as 
iV(A|y| > 0) - N(A\y\ < 0) = f dr f cLy! dy N(r,Y,y)x 

J all J all Jail 

0{y) - 6(-yj\ \e(Y - y)6(Y + y)- 6(-Y + y)0{-Y - y) 



+ 



9(Y) - 6{-Y) 9{-Y + y)6(Y + y) - 9{Y - y)9{-Y - y) 



(5.2) 



We note here that, if we simply define 

A|yf = (y t - yt)9(Y) - (y t - yt)O(-Y), (5.3) 
instead of A|y| = \yt\ — then the numerator of eq.(1.3) becomes, instead of eq.(5.2), 



JV(A|y|' > 0) - N(A\y\' < 0) = f dr f dY ! dy N(r,Y,y)x 

J all J all J all 



e(y)-6(-y) 6(Y)-6(-Y) 



(5.4) 



and we confirm that inclusive values of eq.(5.2) and eq.(5.4) are numerically the same, 
although y- and Y- distributions are different. From eq.(5.4), we find that the numerator 
of the charge asymmetry Ac eq.(1.3) can be obtained by the projection method of eq.(4.6) 
with the weight function of 

w(t,Y) = 9(Y)-9(-Y), (5.5) 

instead of the optimal weight function D A q (r,Y). 

In order to compare the efficiency of the charge asymmetry Ac eq.(1.3) and that of 
the optimal observable eq.(4.6), we introduce a variable N® pt ', 



dr I dY I dy 

< all J all 



6{y)-e{-y) D q j>(T,Y)N(r,Y,y), 



and define the optimal charge asymmetry A r pt ' as 



A°J\i 



(5.6) 



(5.7) 



N T otal(Tj)' 

We evaluate the efficiency of these two observables by using MC generated tt event 
samples assuming the anti-triplet diquark model of eq.(2.1) with = 500 GeV and A = 
2.65. The result is shown in Figure 7. The blue diamonds give the charge asymmetry Ac 
of eq.(1.3) whereas the green squares give the optimal observable A^ 1 ' defined by eq.(5.7). 
The left panel shows the charge asymmetries as functions of the invariant mass of the tt 
system, M#. The right panel shows the distributions of the significance, Asymmetry/Error, 
also as functions of M t j. Only statistical errors are considered. In both panels, each bin 
size of M t i is 100 GeV and the last bin includes all contribution of M t i > 1300 GeV. The 
statistical error of N® pt '(rj) is estimated as 




%) + K-y)] ( D a(t, Y)) 2 N(t, Y, y) . (5. 
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Figure 7. The charge asymmetries (left) and their significances (right) as functions of the tt 
invariant mass, M t j. The blue diamonds give the charge asymmetry of eq.(1.3) whereas the green 
squares give the optimal observable defined by eq.(5.7). Only statistical errors are considered. In 
both figures, each bin size of M t j is lOOGeV and the last bin includes all contribution of M t i > 1300 
GeV. 

Although the significance of the optimal observable A® pt ' is better than that of the 
charge asymmetry Ac in all regions of M t j as expected, they are almost comparable. We 
therefore conclude that the charge asymmetry Ac is a very sensitive observable of the 
sub-process FB asymmetry at the LHC and that it is difficult to require it further once the 
PDF uncertainty are taken into account. 

6 Conclusion 

In this paper, we have studied contributions of a scalar diquark particle in a color repre- 
sentation of anti-triplet and sextet to the top quark production at the Tevatron and the 
LHC. 

We have explored the parameter space of the diquark models which are consistent with 
the measurements both at the Tevatron and at the LHC by including contributions from 
the single and pair production of diquarks. We find that the sextet diquark model has 
a wider allowed parameter space than the anti-triplet diquark model does, and that the 
anti-triplet model can be excluded by the constraint on the charge asymmetry Ac at 90% 
C.L. from the LHC. At 95% C.L., both models still have allowed parameter spaces, which 
can soon be explored at the LHC. 

In addition, we have introduced the optimal observable of the sub-process FB asym- 
metry and compared its statistical significance with that of the charge asymmetry Ac 
measured at the LHC. We find that they are comparable and that Ac is a very sensitive 
measure of the FB asymmetry in the qq — > tt sub-processes. 
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